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We explore some ways of minimally modifying the neutrino mixing matrix from 
tribimaximal, characterized by introducing at most one mixing angle and a CP vio- 



D , lating phase thus extending our earlier work. One minimal modification, motivated 

to some extent by group theoretic considerations, is a simple case with the elements 

^> ■ V a 2 of the second column in the mixing matrix equal to Modifications by 

C\ . 

, keeping one of the columns or one of the rows unchanged from tri-bimaximal mixing 

CO 

all belong to the class of minimal modification. Some of the cases have interesting 
■ experimentally testable consequences. In particular, the T2K and MINOS collabora- 

tions have recently reported indications of a non-zero #13. For the cases we consider, 
the new data sharply constrain the CP violating phase angle 5, with <5 close to (in 
$— 1 ' some cases) and it disfavored. 

PACS numbers: 



I. INTRODUCTION 

Mixing of different neutrino species has been established by various experiments Re- 
cently the T2K Q and MINOS j^J collaborations have reported indications of a non-zero 
#13 providing new evidence for neutrino mixing and new information about mixing pattern. 
The mixing can be represented by the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) 4| mix- 
ing matrix V in the charged current interaction of the W boson with left handed charged 
leptons II and neutrinos vl, C = —{g/\/2)l L ^Vi> L W~ j s [ + h.c. The elements of the unitary 



2 



matrix V are usually indicated by V a j with a = e, /i,r, ... and j = 1,2,3,.... With three 
neutrinos, there are three mixing angles, one Dirac phase, and possibly two Majorana phases 
if neutrinos are Majorana particles. 

It is possible to fit data from various experiments js| , except for possible anomalies in the 
LSND, MiniBoon fl, and MINOS 

a, 

data. The pre-T2K data, with modified Gallium cross- 
section used by SAGE collaboration jg], on mixing and mass parameters can be summarized 
at ler(3er) level as {5, 9|, 

Amii = 7.59 ± 0.20(1^) x 10~ 5 eV 2 , 



Am 2 31 



-2.36 ±0.11(±0.37) [+2.46 ± 0.12(±0.37)] x iO" 3 eV 2 
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Here the mixing parameters are those of the Particle Data Group (PDG) parametrization 



ization [l| 
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Qfj.3, 9 
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9 e $). There is no direct experimental information on the 



(with 9 12 = 9 e2 , 
phases Spdq- 

The mixing pattern is well approximated by the so-called tri-bimaximal mixing pattern 



of the form 
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Note that with this mixing pattern, #13 is equal to zero. 

The recent T2K data on 6* 13 (and therefore (V^l) show that [2j at 90% C.L. it is not 
zero with sin 2 2# 13 in the range of 0.03(0.04) ~ 0.28(0.34) for normal (inverted) neutrino 
mass hierarchy. Data from MINOS ^ also disfavor 6* 13 = at the 89% C.L. level. This 
information provide an important constraint on theoretical model buildings for neutrino 



mixing 



ll| . Combining the recent T2K and M 



NOS data with previous neutrino oscillation 



data and using the new reactor flux in ref.[12j, more stringent constraints on the mixing 
parameters than those given in eq.([T]) have been obtained 13j . These new constraints are 
shown in Table [H We will use them for our later discussions. 

The combined data show that #13 is non-zero at 3a level. In fact the pre-T2K data already 
provides a hint that #13 is non-zero at a more-than-la level [9|. A non-zero V e 3 would rule out 
tri-bimaximal mixing. Now the tri-bimaximal mixing prediction for the angle 9u is outside 
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Parameter 


sin 2 9\2 


sin 2 023 


sin 2 6\3 


Best fit 


0.312 


0.42 


0.025 


lcr range 


0.296 - 0.329 


0.39 - 0.50 


0.018 - 0.032 


2a range 


0.280 - 0.347 


0.36 - 0.60 


0.012 - 0.041 


3cr range 


0.265 - 0.364 


0.34 - 0.64 


0.005 - 0.050 



TABLE I: Ranges for mixing parameters obtained in Ref. 



the la. Also, because one of the elements in the mixing matrix is zero, the tri-bimaximal 
mixing does not allow CP violation to be manifest in neutrino oscillation, i.e. the Jarlskog 



parameter [l^| J is identically equal to zero. 



15 



19] 



will be able to measure CP violation in neutrino oscillation |20|. Were CP violation to be 
established in the future, it would definitively rule out tri-bimaximal mixing. We also recall 
that CP violation in the lepton sector has profound implication regarding why our universe 
is dominated by matter. There is thus an additional motivation to study deviations from 
tri-bimaximal mixing. 



Phenomenologically, small deviations from the tri-bimaxima 
rameterized in terms of three small parameters and studied 21 



pattern can be easily pa- 



22|. In [23|, we studied a 



minimal modification with one complex parameter. Here we extend this discussion. 

Theoretically, many attempts have been made to derive the tri-bimaximal mixing, but 
in our opinion a simple and compelling construction is still sorely lacking. Many of the 
attempts were based on the tetrahedral group first studied by Ma and Rajasekaran 



for neutrino models 
given in 



24J, and subsequently by others [25|. A group theoretic discussion was 
attributing the difficulty to a clash between the residual Zi in one sector and Z3 
in another. Residual discrete symmetries in the context of neutrino mixing have also been 
extensively discussed in Ref. 27J. Authors in Ref. [281] have named this the "sequestering 
problem" . Within the context of A4, it was shown [29| that two assumptions were needed to 
obtain an one-parameter family of mixing matrices which contains tri-bimaximal mixing. In 
other words, to obtain tri-bimaximal mixing, it is necessary to find one reason or another to 
set this particular parameter to zero. This suggests, or at least motivates, studying various 
one-parameter modifications to tri-bimaximal mixing. 
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II. THE FORM OF MINIMAL MODIFICATIONS FOR V TB 



In our 2006 work 



231 ]. which we will review briefly in the appendices, we were naturally 



led in the context of A4, assuming that the sequestering problem of [28[ could be solved, 
to a modification of tri-bimaximal mixing in which the middle "(1,1,1)" column was left 
unchanged. 

It has not escaped the notice of many authors that the 3 columns, 
(—2, 1, 1), (0, 1, —1), (1, 1, 1), of Vtb furnish the 2-dimensional and the 1-dimensional rep 



resentations of the permutation group 5*3. Historically, this observation led Wolfenstein|30] 
long ago to guess, based on a sense that somehow (1, 1, 1) is special, a mixing matrix that 
consists of Vtb with its last two columns interchanged. Another early guess was put forward 
by Yamanaka, Sugawara, and Pakvasa 3l|. The mutual orthogonality of these 3 columns of 
course also mean that they correspond to the three Gell-Mann diagonal generators of SU(3). 



Similarly, these 3 columns also appear 29 



3211 in a table of Clebsch- Gordon coefficients for 



SU(2). These may all indicate some deeper reasons for tri-bimaximal mixing, such as the 



possibility 



29[ that neutrinos are composite. 



With so little known about the underlying theory of neutrino mixing, we take here a 
purely phenomenological approach. If we take the column vectors as reflecting some basic 
feature of neutrino, a minimal modification may be to keep the columns vectors as much as 
possible. Some phenomenological implications have been studies in Ref. 33J. In this paper 
we analyze these minimal modification to the tri-bimaximal mixing in light of the recent 
T2K data. Because of unitarity, we can leave at most one column unchanged. Minimal 



modification to the tri-bimaximal mixing can therefore be characterized by w 



rich co 



umn we 



leave unchanged. Some special cases have been considered in the literature |23l.l29l.l32l.l34j.l35|. 
We refer to them clS CclSG V a , case V\ and case V c } corresponding to keeping the third, second 
and first columns unchanged, respectively. This class of modifications can be obtained by 
multiplying a two generation mixing matrix from the right of Vtb, 2 and 3, 1 and 3, and 
1 and 2 mixing patterns. These modifications can be viewed as perturbation to the tri- 
bimaximal by modifying the neutrino mass matrix 33]. One can also motivate minimal 
modification by perturbing the charged lepton mass matrix in a similar fashion which would 
result in one of the rows ol Vtb unchanged. One such an example has been studied by 
Friedberg and Lee in 35|. Therefore there are another three types of minimal modification, 
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keeping the first row or the second row or the third row indicated by case W a , case W b and 
case W c , respectively. 

A. One of the Columns in Vtb unchanged 



Leaving one of the columns in Vtb unchanged, we have the three possibilities: 



V a = V 



TB 



cost sinr 
— sin r cos r 
1 



V b = Vi 



TB 



cosr sin re* 

1 
smre~ lS cost 



/ 



V c = V 



1 











\ 



V 



cost smre lS ■ (3) 
— sin re~ lS cos r 
We will use the abbreviation c = cos r and s = sin r. 

Note that for V a , with the third column in Vtb unchanged, V e 3 = and there is no 
non-removable phase leading to a vanishing Jarlskog parameter, J = 0. No CP violation 
phenomena can show up in oscillation related processes. 



B. One Of The Rows In Vtb unchanged 



Keeping the neutrino mass matrix unchanged and perturbing the charged lepton mass 
matrix for the tri-bimaximal mass matrices, one obtains a V in the form of a unitary U 
multiplied from left to the tri-bimaximal mixing V = UVtb- The three cases with one of 
the rows unchanged from tri-bimaximal mixing can be written in the following forms 

/ 1 n n \ / « n aoid 



W a 



1 

c s 
-s c 



V 



TB 



-se 



c se 1 
1 
5 c 



V 



TB 
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V 



c 


se ld 





— se~ l 


5 c 











1 



V, 



TB 



(4) 



Since there is no non-removable phase in W a , no CP violation phenomena can show up in 
oscillation related processes for this case. 
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As mentioned earlier, the case W c has been motivated theoretically and studied by Fried- 
berg and Lee [35 ]. 



III. PHENOMENOLOGICAL IMPLICATIONS 



As was mentioned in the introduction, of these six cases, we are theoretically prejudiced 



in favor of V b , which we studied in Ref. 



23] . As far as we know, some of the other cases are 



not well motivated and we analyzed them merely as "strawmen" to be knocked down by 
future precision experiments. 

There are no phases for cases V a and W a . They are given by 



( _2£ 

VE 



ya 
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VE 







V 



w a 
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VE 



1 



(5) 



s c 

' Vs Vz \75 Vs ' V2 V2 

s 

VE Vs V3 VE V? / \ VE T VE Vs V3 V2 ^ Vz / 

The modification to tri-bimaximal is represented by a non-zero s = sinr (c = cost). 
For case V a , the main predictions of this mixing pattern are that |V^3| = l/y/2 which is 
at the boundary of the la allowed range. Present data constrain the modification mixing 
parameter s (with c > 0), which can modify the value for V e 2, to be in the range —0.04 ~ 
0.002(-0.075 ~ 0.0) at la (3a) with the best fit value of -0.02. For case W a , a definitive 
prediction for this case is that V e 2 = l/\/3. This is outside the la range, but within 2a 
range. A non-zero s will modify V^ 3 away from l/y/2. The current data allow s to be in 
the range ~ 0.08(-0.14 ~ 0.105) at the la(3a) with the best fit value to be 0.08. For 
both cases V a and W a , V e % are predicted to be zero. They are in conflict with data at the 
3a level. Also there are no phases for V a and W a . There is no CP violation in neutrino 
oscillation. This provides another test for cases V a and W a . Should future experiments find 
CP violation in neutrino oscillations, these two cases would be ruled out. 

The other four cases have two parameters, one can use available data on the magnitude 
of the elements in V to constrain them and to predict other observables, in particular 
the Jarlskog CP violating parameter J = Jm^iV^V^V^)- Complete determination of 
parameters related to neutrino physics include the mixing angles and the CP violating Dirac 
phase, and also the absolute neutrino masses and possible Majorana phases. Since not much 
information can be used to constrain the Majorana phases, in the following we will use the 
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combined pre-T2K and the recent T2K and MINOS data to study some phenomenological 
implications for the mixing parameters of the other four minimal modifications described 
in the previous section. 



Case V b 



In this case the mixing matrix V b is given by 

/ 

V b = 



2c 
V6 



c 



i 

1 



2se' d 
V6 



V2 V3 V2 VE 



i 



e 

V2 



(6) 



V2 V3 

A definitive prediction of this mixing pattern is that V e 2 = l/y/3. This is outside the la 
range, but consistent with data within 2a level. Since V e2 is fixed to be l/y/3, \V e3 \ can be 
expressed as a function of V e \ with 



V e3 = y/2/3- |Ki| 



(7) 



This can be used as an additional check. 

In this case there is room for a non-zero and also a non-zero Jarlskog parameter given 
by J = | V^i 1 1 V^3 1 sin 5/2 a/3 for CP violation. Using the allowed range for #13, we can easily 
obtain information for s with the best fit value given by 0.177 and the la (3a) allowed range 
0.16 ~ 0.22(0.09 ~ 0.27). 

Expressing V e3 and in terms of the mixing angle r (through c and s) and CP violating 
phase, we have 



e3 



—se 
3 



is 



(8) 



One could use eq. 



to eliminate r and relate \V e3 \, \V„a\ and S 



>3| 



v/3 



cs cos 5 



1/2 



2 1 l / 2 

Q|K 3 | cos5 - Vl - 3|K 3 |72J +^|K 3 | 2 sin 2 5 



(9) 



Using the known ranges of sin #23 and sin #13, the la range for is determined to be in 
the range of 0.617~ 0.701. We could plot |V e 3| in terms of the unknown 5 for some typical 
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values of in the allowed range. In solving for \V e ^\ we have to pick the branch consistent 
with eq.® of course. We will work in the convention where s and c are all positive. The 
result is shown in FigJTJ \V e s\ is symmetric in the region of to tt and it to 2ir for 5. The 
allowed ranges in Table [I] rule out some portion of allowed 5. Regions of 5 close to it are 
not allowed at the la level, but there are ranges of 5 which can be consistent with data. 
Improved data can further narrow down the allowed range. 

Since in this case there are only two unknown parameters, s and S, in the model, the four 
parameters #12, #23, #13 and the CP violating parameter Spdg of the general parametrization 
are related. The boundaries in Fig. [1] represent the la allowed ranges for the parameters of 
the model. The situations are similar for the other 3 cases in our later discussions. Note also 
that the analysis we are carrying out is insensitive to whether the neutrino mass hierarchy 
is normal or inverted. 

0.35 
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0.05 
0.00 

12 3 4 5 6 

s 

FIG. 1: Case V b . \V eZ \ as a function of 5 for |V^3| equals to 0.617 (la lower bound, dashed), 0.641 
(best fit value, solid) and 0.701 (la upper bound, dotted). The solid and dashed horizontal lines 
are for the best value and the la bounds of \V e s\, respectively. 
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One can extract useful information on the CP violating parameter J by using the relation 

J = -i^sin^lKslVl - 3|^ e3 | 2 /2 - (10) 

Note that J is simply proportional to sin 5. The results are shown in Figj2j On the left of 
Fig 12} J is plotted as a function of 5 for three values of (V^l- The absolute value of J can 
be as large as 0.04 for jV^I takes its la upper value of 0.179. 
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FIG. 2: Case V b . 


Figure on the left: J as a 



0.70 - 
0.65 



(dashed), and 0.134 (dotted). Figure on the right: |V^3| as a function of 5 for \V e s\ equals to 0.179 
(solid), 0.145 (dashed), and 0.134 (dotted). The solid and two dashed horizontal lines are for the 
best value and the la bounds of (V^l, respectively. 




FIG. 3: Case V b . Contours of |V^3| and |V^3| for different values of J. The curves are for J equals 
to ±0.04 (solid), ±0.03 (dashed) and ±0.01 (dotted). The solid and two dashed horizontal lines 
are for the best value and the la bounds of |V^ 3 |, respectively. 

For |V^ 3 | close to its lower bound, S close to and 2tt are favored. For 5 close to ir, |V^ 3 | 
is outside of its la allowed range. There are overlaps for the regions allowed in the right 
figure of Fig|2] and those in Fig{TJ Combining information from V e 2 discussed earlier, we can 
conclude that this case is ruled out at la level, but is consistent with data at 2a level. 

In Figj3j the contours of \Ves\ and |V^| for different values of J are shown. The contours 
are degenerate in ±| J|. 



Case V c 
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For this case, the mixing matrix is 



yc 



1_ _c se~ iS 

1 c < se 



se 
V3 



se 1 " 
V3 



V2 V2 
se~^ _ c 



(11) 



A prediction for the mixing is that V e \ = 2/\/6. V e3 is related to V e2 by 



IK 



c3 



Vl/3 -IK2I 



(12) 



Imposing unitarity of the VpmnSi this is prediction in the la region. 

Within the la range of V e3 , s is allowed to vary from 0.2 to 0.31 with the best fit value 
of 0.25, and V e 2 can be within its la range. One can again combine information from \V^\ 
and \V e3 \ to constrain the CP violating phase S and the Jarlskog parameter J. We have 

1 1/2 



K 3 \ 

J = 



V e3 \ cos 5 + -^=\fl 
V2 



3V2 



sin 5 \V f 



3lKs| 2 



3IK 3 | 2 ] | + |K 3 | 2 sin<5 2 



(13) 




FIG. 4: Case V c . \V e3 \ as a function of S for |V^| equals to 0.617 (dashed), 0.641 (solid) and 0.701 
(dotted). The solid and dashed horizontal lines are for the best value and the la bounds of \V e3 \, 
respectively. 

The results are shown in Figs. IU and El For this case even a larger portion of 5 range 
is ruled out by data as can seen from FigJH From the right figure of Fig|3J we also see that 
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FIG. 5: Case V c . Figure on the left: J as a function of 5 for \V e 2,\ equals to 0.179 (solid line), 0.145 
(dashed line), and 0.134 (dotted line) for case V c . Figure on the right: \V^\ as a function of 5 for 
\V e s\ equals to 0.179 (solid), 0.145 (dashed), and 0.134 (dotted). The solid and dashed horizontal 
lines are for the best value and the la bounds of |V^3|, respectively. 
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\vj 

FIG. 6: Case V c . Contours of \V e z\ and |V^3| for different values of J (right). The curves are for 
J equals to ±0.04 (solid line), ±0.03 (dashed line) and ±0.01 (dotted line). The solid and dashed 
horizontal lines are for the best value and the la bounds of |V^3|, respectively. 

a large portion of 5 is ruled out. However, there are overlaps for the two regions. This case 
can be consistent with data at the la level. CP violating information is shown on the left 
in FigJSJ The largest value of J is 0.04 for (V^l equals to its la upper bound of 0.179. The 
correlations of J, V^ 3 and V e3 are shown in Fig® 

Case W b 
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The mixing matrix W b is given by 



/ 



W" 



2c se~~ c I se~" se~~ 

fc fc /o" ' fn ft\ 



1 

2se~ 



V3 y/3 
1 



\ y/6 



V2 
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~V2 ) 



(14) 



This case predicts = l/v2. This prediction is at the boundary of la allowed range. 
But it is different than case V a since V e % is not zero and J can be non-zero, s is in the range 
of 0.19 ~ 0.25 at la level with the best fit value given by 0.20. In this case correlations of 
| V^2 1 , | V^3 1 and 5 or J are given by 

,1/2 



|K 2 | 
J = 



2\V e3 \cos5+ v 7 ! ~ 2|Ks| 2 ) +2|V; 3 | 2 sin 2 <5 



3^ 



sin 5 1 V e 



(:3 



1-2|K 3 | 2 



(15) 



The results are shown in FigsJTJ reffig8 and |9j Now \V e 2\ is playing the role of |V^3| for 
the cases V b and V c . From FigJ7]and the right figure in Fig^, we see that the allowed range 
for 5 is constrained to be even more closer to ir/2 and 3ir/2. But there are still regions 
consistent with data at the slightly larger than la level. J can be as large as 0.04 for \V e s\ 
taking its la upper of 0.179. The correlations of J, and V e s are shown in Figj9] 




1 



FIG. 7: Case W b . \V e 3\ as a function of 5 for \V e 2\ equals to 0.538 (dashed), 0.553 (solid) and 
0.568 (dotted). The solid and dashed horizontal lines are for the best value and la bounds of \V e 3\, 
respectively. 
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FIG. 8: Case W b . Figure on the left: J as a function of 5 for \V e 3\ equals to 0.179 (solid), 0.145 
(dashed), and 0.134 (dotted). Figure on the right: \V e 2\ as a function of 5 for \V e 3\ equals to 0.179 
(solid), 0.145 (dashed), and 0.134 (dotted). The solid and dashed horizontal lines are for the best 
value and the la bounds of (V^l, respectively. 
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FIG. 9: Case W b . Contours of \V e s\ and \V e 2\ for different values of J. The curves are for J equals 
to ±0.04 (solid), ±0.03 (dashed) and ±0.01 (dotted). The solid and dashed horizontal lines are for 
the best value and the la bounds of \V e 2\, respectively. 



Case W c 



In this case the third row is left unchanged from tri-bimaximal mixing 35|and the mixing 
matrix W c is given by 

/ 2c se iS c i se iS se iS \ 
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(16) 
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V2 / 
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We have 



V2\V e3 \cos5 



1 - 2 7, 



c3 



nl/2 



2 IK, 1 2 sin 2 5 



^ = -^sin5|K3lv / l-2|V;3| 2 - 



(17) 



The expression for \V e 2\ is the same as for case W b , but the sign of J is minus of those 
for case W b when reading Figs]8]and|9l respectively. One can easily read off the constraints 
on S and J from FigsJT] and [HJ A crucial difference is that is no longer 1 / y/2, but 
c/V2. \V,a\ can vary from 0.684 ~ 0.694(0.670 ~ 0.703) at the la (3a) level. This 
case is consistent with data at la level. Precise measurement of | V^3 1 can be used to 
distinguish these two cases. 



IV. SUMMARY 



Recent data from T2K and MINOS show evidence of a non-zero V e3 at 90% CL. level. 
There may be the need for modifications to the tri-bimaximal mixing with which V e3 is 
equal to zero resulting in a CP conserving mixing. We have studied several possible ways 
to minimally modify the tri-bimaximal mixing by keeping one of the columns or one of the 
rows in the tri-bimaximal mixing unchanged. Six cases were studied. Two of the cases 
have V e 3 = 0. These two cases are in tension with data at 3a level. Also for these two 
cases, the CP violating Jarlskog parameter J is identically zero. CP violation in neutrino 
oscillation can provide new test for these two cases. For the other four cases, all have two 
parameters in the mixing matrix. Current data on neutrino oscillation can put constraints 
on the parameters, but consistent with data within 2a for case V b , and within la for the 
other 3 cases. The allowed ranges for CP violation are also constrained. Future experiments 
can test the predictions to rule out these models. 
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Appendix A 



Here we review some features of the family group A 4 24| that led us to favor V b , which 



we analyzed in 



23( | . A key point to achieve the tri-bimaximal mixing in models based on 



1 OJ 2 OJ 

1 1 1 

1 OJ oj 2 



\ 



A4 symmetry is to obtain the matrices Ui and U v which diagonalize the charged lepton and 
neutrino mass matrices Mi and M v , UjMiU r = D\ and U^M U U U = D u in the following 
forms 23] 

( \ -1 ^ 
y/2 , (Al) 
10 1 , 

where oj 3 — 1. The mixing matrix V is given by UjU u . Note that U r which plays a role in 
diagonalizing the charged lepton mass does not show up in V. With suitable choices of phase 
conventions, V can be written in the form in eqj2j In general the elements in the diagonal 
matrix D v have phases, the Majorana phases. We provide some details in Appendix B for 
obtaining the relevant mass matrix. 

One can easily find that the neutrino mass matrix must be the following form 



/ 



V 



\ 



7 



( a 0^ 



7 

/3 a 



(A2) 



to obtain the right form for U„. With appropriate redefinition of phases and mixing angles, 
one obtained the form Vtb for V. 

To achieve the above with specific models, there are some requirements for Higgs boson 
fields vacuum expectations values as can be seen in Appendix B. In general, M v above will 
be modified. If one keeps the charged lepton mass matrix unchanged, the most general Higgs 
potential may not respect the conditions leading to modification resulting modifying the '11' 
entry to a — e and the '33' entry to a + e for the neutrino mass matrix M v given by eq JB4l 
in Appendix B. With appropriated redefinition of phases and angles shown in Appendix B, 
one obtains V b where we have redefined 5 = t] + n jl to the form in eq JB8l 

Note that in taking the form of V b in eqj3j without Majorana phases in V, we have 



absorbed possible Majorana phases in the masses fhi 



rrije 



Here m~ are real and 



16 



positive. Without loss of generality, one can always choose one of the «j to be zero, for 
example n 2 = 0. The above form belongs to the minimal modification to the Vtb mixing 
pattern specified earlier with the elements in the second column remain to be V a 2 = 1/ 

We have arrived at this minimal modification from a specific model. In fact, the above 
parametrization is the most general one for V a 2 = up to phase conventions. One 

can understand this by starting with a most general parametrization used by the Particle 
Data Group, Vpdg an d then set certain angles to some particular values to make sure that 
V a 2 = 1/ Vs. Since we want all of the elements in the second column to be 1/ v3, for 
convenience we exchange the second and the third columns, and move the Dirac phase at 
different locations, by redefining the phase of charged leptons and neutrinos, according to 
the following 

V b = P L V PDG EP R , (A3) 

where Pl,r diagonal phase matrices with elements, Pl = diag(e lS , 1,1), Pr = diag(e~ lS , 1,1) 
and E is a matrix switching the second and third columns with elements = 5n5ji + 

We have 



yt 



^ C13C12 S13 Ci3Si2e J<5 

— Sl2 c 23 e lS — C 12 S23 S 13 S23C 13 C12C23 — S^S^iSi^e 15 
\ Sl2S 2 3e~ l<5 - C12C23S13 C23C13 -C12S23 - Si2C23Sl 3 e 1 ' 5 / 



(A4) 



To have all of the elements in the second column in the above matrix to be l/\/3, one just 
needs to set S13 = 1/ \/3, C13 = a/2/3 and S23 = C23 = 1/ V2. We then obtain eq|3] with 
c = C12 and s = su- 

In this basis, one can reconstruct the elements in the neutrino mass matrix in terms of 
the mixing angle r, phase 5 and masses rhi by requiring M v = UfVD u V T llJ . We have 

a= l - ((c 2 - sV M )mi " (c 2 - s 2 e 2l5 )m 3 ) , 

P = \ ((c 2 + #e-"*)m x + (c 2 + s 2 e 2lS )m 3 ) , 
e = ics(e~ tS fhi — e zS m 3 ) , 7 = m 2 - 

The a, (3 and e here are equivalent to those given in Appendix B with a different basis. The 
important thing is that the number of independent parameters are the same, total six of 
them. 
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Appendix B 



In Appendix B, we briefly outline how V b can be obtained from models based on A4 



family symmetry 24|] with the standard model (SM) gauge symmetry. 

A way to obtain XJ\ in eq JAll is to assign the three left handed SM lepton doublet II = 
(Ilx, Il2, Il3), and three right handed neutrino {Iri, Ir2, Irs) representations into a triplet, 
and the three singlets 1, V and 1" of A 4 group, respectively A SM doublet and A 4 triplet 
Higgs representation $ = ($1, $2, $3) then leads to the Yukawa coupling terms 



U = \ e (lL$i)lm + A M w 2 (/ L $)i// R2 + \ t u(Il®)i»Ir3 + h.c. 



(Bl) 



After $ develops a vev of the form < $ >= (v<f>,v$,v<f>), the charged lepton mass matrix is 
given by, Mi = U1D1 with D\ = Diag(m e , m^, m T ) and m ; = \/3v§\i. This gives the right 
U\ with U r = I. 

To obtain the right form of neutrino mass matrix, three right handed neutrinos Vr = 
{vri, vr2, vrz) of A4 triplet, a SM doublet and A4 singlet Higgs 0, and a SM singlet and 
triplet Higgs \ — (Xi? X2, Xs)i are needed. The Yukawa coupling terms relevant are 



L v = K(I l vr)i + m^RV^x + \ x {{vrVr)zX)i + h.c. 



C\ 



(B2) 



With the vev structure < <ft >= and < \ > = (0, v x , 0), the neutrino mass matrix is given 
by 





f 


Md \ 




M= | 




, Mr = 




\M T D 


Mr J 



( 



m 
m 



y m x m 



\ 



J 



(B3) 



where Mo = Diag(l, 1, l)X v v,p, and m x = X x v x . The light neutrino mass matrix is M v = 
—MdM^M^ which is the desired form giving U u in eq iAll 

The charged lepton masses are controlled by the vev of $ which requires the components of 
<3> to have the same vev which leave Z% unbroken symmetry in the theory, while the neutrino 
mass matrix is controlled by the vev of <fi and x with only < X2 > non-zero preserving a Z2 
unbroken symmetry. If there is no communication between these two sectors, the residual Z 3 
and Z2 symmetries are left unbroken. But in general these two sectors cannot be completely 



sequestered and can interact which complicates the situation |26l. 



28 



for example the 
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vev structure with only < X2 > De non-zero may not be maintained. One of the 

consequences concern us is that the mass matrix M v will be modified to [23] 



( a-e ^ 



M v 



\ 



7 
P a + e 



(B4) 



The above will lead to a different form for U v from that in eq JAll which can be written as 



where cq = cos 9 and sq = sin 9 with 

4|/3| 2 



-cq -is e ^ 
1 
—see tp icee jp 



(B5) 



tan 2 (2#) 



1 



Ala 



2 £ 2| 



(|a + e| - |a - e|) 2 (\a + e| + \a - e\) 2 

tan a 



sin 2 cr) , <7 = arg( 



P 



' OL — 6 

p = 5 + arg(<y/ — ; — ), tan 5 



a + e 



a - 


-e| - 


- \a 


— e 


a - 


he| - 


- \a 


— e 



(B6) 



The phase conventions are chosen such that setting e = 0, the resulting mixing matrix V 
goes to Vtb- 

In this basis, V is given by 



r 



i 



v/3 



— (wee + u 2 see lp ) 1 i(u 2 cge ip — usg) 

— (u 2 cg + ws6ie 4p ) 1 i(ucge lp — oj 2 se) 



(B7) 



which can be further rewritten as 



V = Vn 



TB 



( 



cos r i sin re* 77 
1 
i sin re~ ir? cos r 



(B8) 



where V p is a diagonal phase matrix V^, = (e*^ +p / 2 \ l,e^ £+^/ 2 ) multiplied from right, r, 77 
and £ given by 



sin 2 r = -(1 — sin(20) cos p) , tan£ 



1 — tan# 
1 + tan 9 



tan(p/2) , tan?7 = tan(2#) sin p . (B9) 
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Absorbing the phases in V p into the neutrino masses, the total Majorana phases for rhi, 
m 2 and m 3 are k x = — (2£ + p + 2ai), k 2 = a 2 and k 3 = — (— 2£ + p + 2a 3 ) with 

ai = -i[arg(cg|a - e| + 2s e c e \l3\e^ +r7) + s 2 e \a + e\e 2i ' s )] - arg(- v / ^^) , 
«2 = --axg(j) , 

a 3 = --[arg{<%\a - e\ - 2s e c e \l3\e^ +a) + s 2 e \a + e\e 2l ~ 5 )} - arg(z'V^Ti) . (BIO) 

In this basis, V is given by eqfB8j but with V p removed. 
The absolute masses squared are given by 

m ? = |a| 2 + | e | 2 + |fl| 2 2R < ae *) 



cos(2^) 



m 2 = \j\ , 

^=w + M 2 +w 2 +Wr- (B11) 



[1] Particle Data Group, C. Amsler, et al., Phys. Lett. B667, 1 (2008), and 2010 edition. 

[2] The T2K Collaboration, K. Abe et al., larXiv: 1106.28221 

[3] P. Adamson et al., The MINOS Collaboration. larXiv:1108.00151 

[4] Z. Maki, M. Nakagawa, and S. Sakata, Prog. Theor. Phys. 28, 870 (1962); B. Pontecorvo, 

Sov. Phys. JETP 26 (1968) 984 [Zh. Eksp. Teor. Fiz. 53 (1967) 1717]. 
[5] M. Gonzales-Garcia, M. Maltoni and J. Salvado, larXiv: 1001 .45241 

[6] A. Aguilar et al., LSND collaboration, Phys. Rev. D64, 112007(2001); A. A. Aguilar-Arevalo 

et al., MiniBooNE collaboration. larXiv:1007.1150[ hep-ex]. 
[7] P. Vahle, MINOS collaboration, talk at the 24th International Conference on Neutrino Physics 

and Astrophysics (Neutrino 2010), June 14, 2010 (Athens, Greece). 
[8] SAGE collaboration, J.N. Abdurashitov et al, Phys. Rev. C80, 015807(2009). 
[9] G. Fogli et al, J. Phys. Con. Ser, 203, 012103(2010). 
[10] P.F. Harrison, D.H. Perkins and W.G. Scott, Phys. Lett. B530, 167(2002); Z.-Z. Xing, Phys. 

Lett. 533, 85(2002); Xiao-Gang He and A. Zee, Phys. Lett. B 560, 87(2003). 
[11] Z.-z. Xing. larXiv:1106.3244t Nan Li and B.-Q. Ma. larXiv:1106.3284; E. Ma and D. Wegman, 

larXiv: 1106 .4269 1 Y.-J Zhang, B.-Q. Ma, larXiv: 1106.40401 



20 



[12; 

[13 
[14 
[15 
[16 
[17 
[18 
[19 
[20 
[21 
[22 



T.A. Muller et al, Phys. Rev. C 83, 054615 (2011) [arXiyJ10L2663 [hep-ex]] . 

G.L. Fogli et al.. larXiv:1106.60281 

C. Jarlskog, Phys. Rev. Lett. 55, 1039(1985). 



F. Ardellier et al., hep-ex/0405032 



S.-. B Kim (RENO), AIP Conf. Proc. 981, 205(2008). 
Y. Itow et al, Nucl. Phys. Proc. Suppl. Ill, 146(2001). 



X. Guo et al, (Daya-Bay), hep-ex/0701029 



I. Ambats et al. (Noi/A), hep-ex/0503053 



P. Huber, M. Lindner, T. Schwetz and W. Winter, JHEP 0911:044(2009). 
A. Zee, Phys.Rev. D68, 093002(2003) [arXiv:hep-ph/0307323| 

N. Li and B. Q. Ma, Phys. Rev. D 71, 017302 (2005); S. Pakvasa, W. Rodejohann and 
T. J. Weiler, Phys. Rev. Lett. 100, 111801 (2008); S.F. King, Phys. Lett. B659, 244(2008); 
X. G. He, S. W. Li and B. Q. Ma, Phys. Rev. D 78, 111301 (2008); Phys. Rev. D 79, 073001 
(2009); G. K. Leontaris and N. D. Vlachos, !larXiv:1103.6178l [hep-ph]; J.-M. Chen, B. Wang 
and X.-Q. Li, larXiv:1106.3l"33l 



[23] Xiao-Gang He and A. Zee, Phys. Lett. B645, 427(2007), [hep-ph/0607163, (2006). 
[24] E. Ma and G. Rajasekaran, Phys. Rev. D64, 113012 (2001); E. Ma, Mod. Phys. Lett. A 17, 
627 (2002). 

[25] K. S. Babu, E. Ma and J. W. F. Valle, Phys. Lett. B552, 207(2002); G. Altarelli and F. 



Feruglio, Nucl. Phys. B720, 64(2005); K. S. Babu and X.-G. He, |hep-ph/050721~7] E. Ma, 
Phys. Rev. D72, 037301(2005); Phys. Rev. D73, 057304(2006). 
[26] A. Zee, Phys. Lett. B630, 58(2005). 

[27] C.S. Lam. larXiv : 1003.0498 and earlier work cited therein. 

X.-G. He, Y.-Y. Keum and R. Volkas, JHEP 0604, 039(2006). 
P. Kovtun and A. Zee, Phys.Lett. B640, 37(2006). 
[30] L. Wolfenstein, Phys. Rev. D18, 958(1978). 

[31] Y. Yamanaka, H. Sugawara and S. Pakvasa, Phys. Rev. D 25, 1895 (1982) [Erratum-ibid. D 
29, 2135 (1984)]. 

[32] J. D. Bjorken, P. F. Harrison and W. G. Scott, Phys. Rev. D 74, 073012 (2006). 

[33] C. H. Albright and W. Rodejohann, Eur. Phys. J. C 62, 599 (2009) |arXiv:0812. 04361 [hep- 



21 



ph]]; C. H. Albright, A. Dueck, W. Rodejohann, Eur. Phys. J. C70, 1099-1110 (2010). 
|arXiv:1004.2798l [hep-ph]]. 
[34] I. Stancu and D. V. Ahluwalia, Phys. Lett. B 460, 431 (1999) [arXiv:hep-ph/9903408] ; 
C. S. Lam, Phys. Lett. B 507, 214 (2001); C.S. Lam, Phys. Lett. B507, 214(2002); 
W. Grimus and L. Lavoura, Phys. Lett. B 572, 189 (2003); R. Friedberg and T. D. Lee, 



arXiv:hep-ph/0606071 Z. z. Xing, H. Zhang and S. Zhou, Phys. Lett. B 641, 189 (2006); 



R. N. Mohapatra and H. B. Yu, Phys. Lett. B 644, 346 (2007). 
[35] R. Friedberg and T.D. Lee. larXiv: 1008.0453 [ hep-ph]. 



